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1 Introduction 


There are two kinds of numerical methods for option pricing; one is based on the proba¬ 
bilistic approach and another one is the finite difference method for PDE. 

The binomial tree method (BTM), first proposed by Cox, Ross and Rubinstein © , is 
one of the probabilistic numerical methods for pricing options. Due to its simplicity and 
flexibility, it has become one of the most popular approaches to pricing options. mm 


14 18 19 


It is well known that the BTM for European option in Black - Scholes diffusion model 
converges to the corresponding continuous time model of Black and Scholes ( 1) . In 


particular, Jiang 10 showed that the BTM for European option is equivalent to a special 


explicit finite difference scheme for Black-Scholes PDE and proved its convergence using 
PDE approach. 

Amin and Khanna [2] first proved the convergence of BTM for American options using 
probabilistic approach. 

Jiang and Dai ( |l 1,12]) proved the convergence of explicit difference scheme and BTM 
for American options using viscosity solution theory of PDE. They showed that the BTM for 
American option is equivalent to a special explicit finite difference scheme for a variational 
inequality related to Black-Scholes PDE, proved monotonic property of the price by BTM 
and explicit finite difference scheme, existence and monotones of approximated optimal 
exercise boundary and used the method of Barles et al [3][4] and comparison principle 
of |3,7 . Jiang and Dai 13 studied the convergence of BTM for European and American 
path dependent options by PDE approach. Liang et al 16 obtained a convergence rate 
of the BTM for American put options with penalty method and Hu et al 9| obtained an 
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optimal convergence rate for an explicit finite difference scheme and BTM for a variational 
inequality problem of American put options. 

BTM is extended to the jump-diffusion models for option pricing. Amin 1| generalized 
their algorithm of j2| to jump-diffusion models. Zhang 25 studied numerical analysis for 
American option in jump-diffusion models. Xu et al 24 studied numerical analysis for BTM 
for European options in Amin’s jump-diffusion models and gave strict error estimation for 
explicit difference scheme and optimal error estimation for BTM. Qian et al 23 proved 


equivalence of BTM and explicit difference scheme for American option in jump-diffusion 
models, convergence of explicit difference scheme, existence and monotones of optimal ex¬ 
ercise boundary. Luo 19 studied approximated optimal exercise boundary of American 
option in jump-diffusion model. Liang 15 obtained a convergence rate of the BTM for 


American put options in jump-diffusion models. Liang et al 17 obtained an optimal con¬ 
vergence rate for BTM for a variational inequality problem of American put options in 
jump-diffusion models and a convergence rate estimate of approximated optimal exercise 
boundary to the actual free boundary. 

The above all results are obtained under the assumption that the interest rate and 
volatility are all constants. 

On the other hand, Jiang 


10 studied Black-Scholes PDE with time dependent coef¬ 


ficients as a model for European options in diffusion model and provided the generalized 
Black-Scholes formula. H.C. O et al 22 derived a pricing formula of higher order binary 


with time dependent coefficients and using it, studied the pricing problem of corporate zero 
coupon bonds. Such higher order binaries with time dependent coefficients are arising in 
the pricing problem of corporate bonds with discrete coupon ( [21]). H.C. O et al 20 


studied some general properties of solutions to inhomogeneous Black-Scholes PDEs with 
discontinuous maturity payoffs and time dependent coefficients. 

This article concerns with binomial tree methods and monotonic properties for Amer¬ 
ican put options with time dependent coefficients. We consider monotonic properties and 
convergences of prices by binomial tree methods and explicit difference schemes for the 
variational inequality model of American put options with time dependent coefficients and 
then using them prove the decreasing property of the price of American put options and 
increasing property of the optimal exercise boundary on time variable. 

When the coefficients are time dependent, in particular, in the case with time dependent 
volatility, it is not reasonable to assume that the dynamics of the underlying assets price 
forms a binomial tree if we use a partition of time interval with equal parts. Thus one of 
our main problems is to find a time interval partition method that allows binomial tree 
dynamics of the underlying assets price. Another point is to prove the monotonic property 
of option price and approximated optimal exercise boundary. Jiang and Dais convergence 
proof ( [12] ) strongly depends on the monotonic property of option price but such monotonic 
property of option price may not hold when coefficients including interest rate and volatility 
are time dependent as you can see in the following remark 4. We found a special time 
interval partition method and conditions under which the prices of American put option by 
BTM and explicit difference scheme have the monotonic property on time variable. Such a 
special partition of time interval needs some annoying consideration in proving convergence 
to viscosity solutions. 

The remainder of this article is organized as follows. In section 2, we find a time interval 
partition method that allows binomial tree dynamics of the underlying assets price and 
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briefly mention BTM for European options. In section 3, we study BTM price of American 
put option, its monotonic property and existence of approximate optimal exercise boundary. 
In section 4, we study explicit difference scheme for variational inequality model for Amer¬ 
ican put option and show the monotonicity of option price on time-variable and existence 
of approximated optimal exercise boundary. Section 5 is devoted to the convergence proof 
of the explicit difference scheme and BTM. 

2 Time Interval Partition and BTM for European Options with 
Time dependent coefficients. 

Let r(t), q(t) and a(t) be the interest rate, the dividend rate and the volatility of the 
underlying asset of option, respectively. Let 0 = to < ti <••• < tw = T be a partition of 
life time interval [0, T\ and denote as follows 

r n = r(t n ), q n = q(t n ), cr n = cr(t n ), 

Pn = 1 + qn At n ■ Pn = 1 T Af n , 

At n — tra+1 tni — 1, * ' ^ N 1. 

The volatility o(t) of the underlying asset determines the fluctuation of its price in time 
interval [t, t + At]. So if we divide [0, T] by equal parts, then the dynamics of the underlying 
assets price in subinterval [t n ,t n + 1 ] of time may not form a binomial tree. It makes BTM 
difficult in the case with time dependent coefficients . 

On the other hand, from the practical meaning of the volatility <j„, although we consider 
the underlying assets price S in the some interval [t n ,t n+ 1 ] , the underlying assets price S 
largely changes if o n is large; the underlying assets price S changes a little if a n is small. So 
we can imagine that we can make the widths of changes of S in all subintervals a constant 
if we differently define the length A t n = t n+ 1 — t n of subinterval [t„, t n+ 1 ] according to the 
size of cr n . In other words, if we define A t n = t n+ 1 — t n such that a n ■ At n = const = (In u) 2 
, then we can assume that the width of change of S in every subinterval [t n ,t n +i\ is u and 
the dynamics of S in every subinterval 1 ] satisfies one period - two states model 10 . 

Then St are random variables and the evolution in [0, T] forms a binomial tree. Such a 
partition method provides a key to overcome the difficulty arising in the case with time 
dependent coefficients. 

Let us define t n (n = 1, • • • , N ) more definitely. Let assume u > 1. First, we define 

t 0 = 0, <j 0 = cr(f 0 ), Af 0 = - 2 > t 1 =t 0 + At 0 = (lnu) 2 • —. 

a o °o 

If t\ < T, then we define as follows: 

(In it) 2 . , 2 / 1 1 \ 

dr = cr(fi), Ati = - 2 j t 2 = ti + Ati = (lnu) 2 • ^ . 

crj V CT o °i) 

Inductively, if t n < T, then we define as follows: 

. . . (lnu) 2 . ,, , 9 / 1 1 \ 

<y n = o(t n ), At n = - t n+ 1 =t n + At n = (lnu)" • —2 H-b -2 • ( 1 ) 
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Such a process is continued until tjy < T < tjv+ 1 - Then the number N of subintervals 
depends on u,T and a(t). 

If we assume that 

0 < a < a(t) < d, (2) 


then we obtain lower and upper bounds for the size A t n of subintervals of time and the 
number N of subintervals. From the definition (JT|) (2.1) of A t n , we have 


<ln “> 2 < At , < (1 ”" )2 


- 


- /t2 


(3) 


On the other hand, if we use tjy <T < t^+i, then we have 

T-c r 2 i Ar T ■ d 2 

-=-1 < N < -. 

(ln u) 2 - (lnzt) 2 


(4) 


Remark 1. If u J, 1, then N —> oo and 0 < T — tw < A = (lnzt) 2 - A- > 0 . 

Now we consider the dynamics of the underlying asset’s price S. Assume that the width 
of change of S in every subinterval [t n ,t n +i] is u , d = u” 1 and the dynamics of S in every 
subinterval [t n ,t n +i] satisfies one period - two states model. That is, the underlying assets 
price S trl at time t n is changed into S tn u or S tn d . If the initial price of S is Sq , then S tn 
can take one of the following values 

S™ = Sou n ~ a d a (0 < a < n) or Sj = Squ j (j = n,n — 2, • • • , — n + 2, — n). 


Assume that 


dljn < Pn < UTJn , U = 0, 1, • • • , N. 


If we denote 


Pn / Pn - d 

U — d ’ 


n = 0,1, • • • , N. 


(5) 

( 6 ) 


then we have 0 < 9 n < 1 and BTM price of European option with time dependent coeffi¬ 
cients is provided as follows: 


V? = ( S% - E) + (for call) or {E - S*) + (for put), 0 < a < N, 

V£ = —[0 n V2£ + (1 - e n )V^}, 0 < a < n, n = N - 1, • ■ ■ , 1, 0. (7) 

T]n 

Remark 2. Using Jiang’s method ( 101), we can easily prove the followings: BTM can 
be seen as a special explicit difference scheme for Black-Scholes PDE 


dV 

~dt 


a 2 (t) 8 2 V 
2 dx 2 ~*~ 


r{t) - q(t) 


a 2 {t) 1 dV 
2 dx 


r(t)V = 0, oo < x < oo, 0 < t < T, 


V(x,T) = (e x — E) + or (E — e x ) + , — oo < x < oo. 


( 8 ) 











BTM and EDS for American put options with Time Dependent Coefficients 


5 


Let x m = to Ax (—00 < m < 00), 0 = to < fi < • • • < fjv = T and A t n = t n +1 — t n . Denote 
Vff = V(x m , t n ) . Then the explicit difference scheme for ([8]) is provided as follows: 


N 


(e mAl -E)+ or (E - e m * x ) 


mAa;\ + 


1 


1 + r n At r 


1- 


Ax 2 


V n+1 + 

v m 


A At r, 

2Ax 2 


Vn At r, 

2Ax 2 



Atn 
Ax 


yn+l 
m— 1 



Atn 

Ax 


Ktl 


n = N- 1, - - - , 1, 0. (9) 


The scheme © is consistent if r(t) ,q{t) and a(t) are bounded and continuous on [0,T]. 
Such an explicit difference scheme is stable if 


a^Atn < Ax 2 ; 


i-4 


T n Qn 


Ax > 0, 


0 < Vn< N- 1. 


Let Ax —> 0 . Then A t n —> 0 and E converges to the solution to ([s]) . So BTM price 0 
also converges to the solution to 0- 


3 BTM for American Put Options with Time Dependent Coeffi¬ 
cients. 

Let 0 = to<ti<---<t 7 v<Tbe the partition of time defined in (|T]) and let 

Sj = Squ 1 (j = n, n — 2, • • • , — n + 2, — n ; n = 0, • • • , N), 


Vi = ( E ~ s i) + - 

Then BTM prices VJ 1 = V(Sj,t n ) of American put option are provided as follows: 

Vf , = <Pj, 

v.J 1 = max ( — [OnVf+i + (1 - 9 n ) P/.+ 1 ] , Vj \ , n = N — 1, • • • , 1, 0. (10) 

{Pn J 

Now we consider the monotonic property of BTM price VJ 1 for American put option. 


Theorem 1 BTM prices of American put option 

VJ 1 = P(Sj,t n ; E) (n = 0, 1 ,N ,j = n, n - 2, • ■ ■ , —n + 2, -n) 

are decreasing with respect to Sj and increasing with respect to E. That is, 

VJ 1 = P(Sj,t n ; E) > P(S j+1 ,t n ; E) = V? +1 , 

P{Sj,t n ; £d) < P(Sj,t n ; E 2 ) if E x < E 2 . 


( 11 ) 
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Proof Vj* = < pj = (E — Sj) + is decreasing function on Sj and increasing on E. Now 
assume that V k+1 > V^ 1 when n = k + 1 . Then we have 

V? = max|^[0 fc F/ + Y + (l-0 fc )P/_ + 1 1 ],^J 

> max | \e k v^ + (1 - 9 k )V k+1 ] , Vj+1 1 = F/ +1 . 

Thus V k is decreasing on Sj. Similarly, we can prove Vf is increasing on E. (QED) 

In order to prove that V™ is decreasing on time variable, we need the following lemmas. 


Lemma 1 (i) If r{t)/<j 2 {t) is increasing on t, then p n < p n +i- 

(ii) If q(t) / a 2 (t) is decreasing on t, then rj n > rj n+ i. 

(iii) If r[t)/a 2 (t) is increasing, q{t)/a 2 (t) decreasing and A t n is sufficiently small, then 

Pn/'hn A Pn+l/Pn+1 ■ 9 n A 9 n j. 1* 


Proof (i) If r(t)/a 2 (t ) is increasing, then from the definition of At n , we have 


r „+1 


’n+l 


> ^ O (In u) 


2 r n +1 


> (lnw ) 2 




'n+1 


Pn+ 1 — 1 T r n+ iAt n+ i > 1 T r n At n — p r 


(ii) is proved in similar way with (i). 

(iii) If A t n is sufficiently small, then p n > 0 . Since p n < p n + 1 and rj n > r] n+ 1 , we have 
Pn/Vn < Pn+i/Vn+1■ Thus from ([6|, we have 9 n < 9 n+1 . (QED) 

Lemma 2 (i ) If A < B and 0 < a < /3, then aA + (1 — a)B > (3A + (1 — fd)B. 

Proof aA + (l- a)B - /L4 - (1 - /3)B = (0 - a){B - A) > 0. (QED) 


Theorem 2 Assume that © is satisfies, r(t)/a 2 [t) is increasing and q(t)/a 2 (t ) decreasing 
on t. Then for BTM prices V™ of American put option we have 

yn-l > VJ 1 . 


Proof 


From (101 we have 


V J N ~ 1 >ip j = Vf f (j = N,N -2, - ■ ■ , —N + 2, —TV). 
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Now assume that Vj* > V k+1 (Vj). Then we have 

V*~ x = max{— [9k-iVf +1 + (l-0 k - 1 )V!'- 1 \ ,^\ 

l Pk—l J 

> max ( — [0 fc _ iVf+j 1 + (1 - e k . x )V k + x ] , X 

L Pk —1 J 

> max | i- [6» fc _i^*+ 1 1 + (1 - 0 fc _i)V/+! 1 ] , ^| 

> max | [6 k V k + x + (1 - e k )v^ x ] , ^ +1 1 = V* 

Here the first inequality comes from the induction assumption V!f > V^ +1 (Vj) , the second 
inequality from lemma |T] (i), the last inequality from lemma [l] (iii), theorem [l] and lemma 
i (QED) 

Remark 3. Theorem [2] strongly represents the effect of time dependent coefficients. 
Here the main tools are lemma |T] and lemma [2] 

Remark 4. The conditions of theorem [2] are essential. See the following figures: 

V 



Figure 1: Plot (t n ■ V(Sj,t n )) when r(t) = 0.1, q = 0, a = 1, T = 5, E = l,j = 1 


Remark 5. Only using the analogs of lemma |T| and lemma [2] it seems difficult to prove 
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Plot(t n ,Sj,V) when r[t]=0.1,0st<5,q=0,a=l; 

T=5,N=10,At=T/N;u=Exp[cr Sqrt[At] ] 


Figure 2: r(t) is increasing, so V is decreasing on t. 


V 



Figure 3: Plot (t n : V(Sj,t n )) when r(t) = Piecewise{{ 0.2,0 < t < 2}, {0.1,2 < t < 
5 }}i q = 0, a = 1,T = 5, E = 1, j = 1 
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Plot (t n , S 3 , V) when r [t] = Piecewise [{(0.2, 0 s t < 2}, (0.1, 2its5}}], 
q = 0, CT = 1; u = Exp [a Sqrt [At] ] T = 5 , N = 10, At = T / N; 

Figure 4: r(f) is not increasing, so V is not decreasing on t 

that American call option’s BTM price is decreasing on t. 

Now we consider the existence of approximated optimal exercise boundary. 

Theorem 3 Let A t n be sufficiently small. Under the conditions of theorem [1] , for every 
t n {0<n<N~l), there exists a j n £ Z such that 

VJ 1 = <Pj for j < j n , 

VJ 1 > Vj for j = j n + 1, 

VJ 1 < ipj for j > j n + 2. (12) 

Furthermore we have 

jn— 1 — jn- (1^) 

Proof Without loss of generality, we assume that S 0 = 1 and E = 1. (Otherwise, use 
change of variables S = S/E , h = V/E.) Since 

Vf T = (l-S j ) + ='p j , (j = N, N — 2, ■■ ■ , —N + 2, —N), 

we have 

<Pi = Vf f = Q (j > 0); = Vf > 0, (j < -1). 

Since V] N ~ l = max [d N _ 1 (p j+1 + (1 - d N _ 1 )ip j _ 1 ], |, we have 

-1 - o = vj* = Vi, j > 0. 
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In particular, 


V? 


— <Pj (j l)i — Pjv—i(l @N—1)^—1 > 0 — y?o- 


(14) 


Now we consider the case of j < — 1. 


V'j N 1 = max 


= max 


{ —— [9n-i<Pj+i + (1 - Ou-i)<Pj-i] j <Pj i 
L Pn-i j 

— [0JV-l(l - U j+1 ) + (1 - 0JV—l)(l - M- 7-1 )] ,<Pj\ 
lPN-1 J 

= max {pj f 1 _ 1 - Tin-iV 3 , l-u 3 } . (15) 

Note that if j — > —00 , then pj l 1 _ 1 — pj f _ 1 'u 3 < 1 and 1 — u 3 —► 1 . So there exists 
3n- 1 = max{j e Z : j < -1, - Vn-i^ <l-u 3 }. 

If j < jjv-i , then we have pj/_ 1 — V~n_i ^ 1 — u 3 and thus V)' v ~ 1 = 1 — u? = If 

jlv-i + 1 < j < — 1 , then we have pj/_ 1 — Vn-i > 1 — u 3 and thus V^ -1 > ^ . So jJv-i 
satisfies (12) with n = N — 1 . (In particular if ?pv_i < l(<t=> qw-i < 0) , then j N _ 1 = —1.) 

Now assume that when n = k , there exists j k satisfying (12) and (13). Then if j < j k — l, 
then V J k _ 1 = <Pj-i, V k +1 = <p>j+i and thus from the formula (10) and the same calculation 
in (15) we have 


V?- 1 = max 


Pk-l 


[9 k _ ip j+1 + (1 - d k _ 1 )(f j _ 1 ] ,ipj\= max 


1 


u J 


Pk-l Pk-l 


In the case that p k -\ > 1 (•<=> q k -1 > 0) , let 

l = max{j £ Z : j < j k - 1, p^\ - pjj^u 3 < 1 -u 3 }, 


1 — u 3 


then we have l < j k — 1 . If l < jk — 1 j then we define jk- 1 = l ■ Then using the similar 
way with the consideration when n = N — 1 and theorem [2j we have 

j < jk- 1 => Vf -1 = <Pj\ j = jk- 1 + 1 => > V’j! J > jfc-i + 2 => >V k > <Pj. 

If l = jk ~ 1 (that is, V^ -1 = for all j < j k — 1 ), then note that 

j = j fc _! + 1 =► Vf- 1 > Vf > <pf, j > j k - 1 + 2 =► I//" 1 > V* > Vj . 

Generally, we have V k > ipj when j = j k . So if V k k > tpj k , then we define j k -i = jk — 1 ■ 
If Vji = 1 then we define j k -\ = j k . Thus in any case j k - 1 (< Jfc) is well defined. (QED) 
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4 Explicit Difference Scheme for Variational Inequality Model of 
American Options with Time Dependent Coefficients. 

A Variational Inequality pricing model of American option with time dependent coefficients 
is provided as follows: 

min< ! - ( r 0) - q W s % + r(t)V,V = 0, 


0 < t <T, 0 < S < oo, 
V(S,T)=ip{S), 0 < S < oo. 


(16) 


Here 


^(S) = (S-E)+ (for call), = (E - S) + (for put). 

Using the transformation 

u(x,t) = V(S,t)-, S = e x , 

the problem (16) is changed to the following problem 
du a(t) 2 d 2 u 


(17) 


mm <-— 

S dt 2 dx 2 


- r(t) - q(t) - 


cr(t) 2 \ du 


-+r{t)u,u-ip \ =0, 


0 < t < T, — oo < x < oo, 
u(x, T) = ip(x), — oo < x < oo. 


(18) 


Here 


(p(x) = ( e x — E) + (for call), ip(x) = (E — e x y (for put). 


We construct a lattice on E = {—oo <:r<oo,0<t<T} as follows: Select any c £ R 
and Aa; . Let Xj = jAx + c . When 0 < a < 1 , we define as follows: 

rv/\qr 2 rvAx 2 

to = 0, At 0 = ^7—, ti = t 0 + Ato, At 1 = 


’(to)’ 


2 (U)’ 


rv A 

tn — tn — 1 H - 1? — 9 ( ± “5 — 0? 2, 


2 (t n y 


(19) 


This process is continued until tw such that 


aAx 2 

tN = tjV -1 + AtN-1 < 1 < tp {+1 = fjV H- 777 — 7- 

Then we have a lattice on E = {— oo < x < oo, 0 < t < T}: 

Q c = {( Xj,t n ) : Xj = jAx + c, 0 < n < N,j G Z}. 

Under the assumption ([2]) we have 

aAa : 2 . aAx 2 
—< At n < — 


( 20 ) 


( 21 ) 
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Thus there exists TV such that t/v < T < tN+i and we have 


T -a 2 , Ar T ■ a 2 

- —1<N< 


aAx 2 


aAx 2 


( 22 ) 


0 . 


Therefore if Ax —> 0, then TV — > oo and 0 < T — £jv < At jv < 

Uj = u(jAx + c, t n ) represents the value of approximation at (jAx + c, t n ) and let <pj = 
<p(jAx + c). Taking explicit difference for time and the conventional difference discretization 
for space variable in ( 181 , we have 


u] +1 - u] a 2 (t n ) - 2 u] +1 + 


n+1 


A t n 

r{t n ) - q(t n ) - 


2 

& 2 (tn) 


n+1 _j_ 0< n+l 
J _ 

Ax 2 

n+1 n+1 


u j+i ~ u j-i 

2Ax 


+ r(t n )u ”, u” - ifj > = 0. (23) 


If we denote r n = r{t n ), q n = q{t n ), &n = cr(t n ) , then (231 is equivalent to 


Uj = max 


1 -b v n At n 


1 Ax 

+ 2 “ 2 a 2 


1 - 


ra +1 , 

Ax 2 / J Ax 2 


1 Ax 

2 + 2a 2 


r n ~ Qn - 


,n+l 
l j +1 


r n - 


'b-i 




Here, if we denote £>0 = e c , then 

= (Sbe- 7 ^ — E) + (for call), ipj = (E — S'oe :,Ax ) + (for put). 


From (19) we have a = 


Ax 2 


and let 


1 Ax 


2 2a 2 

Then we have the explicit difference scheme 
V? = <Pj, J € z, 


r n ~ Qn- 


Uj = max { — {(1 - a)t /; +1 + a [c a n U. 7+ 1 + (1 - a n )U^] } , <pj ) 

L Pn J 


(24) 


(25) 

(26) 


(Note that p n = 1 + r n At n .) In particular, if a = 


a^ i At n 

Ax 2 


n = N- 1,... , 1, 0. 

= 1 , then Ax = o n \J At n and 


Uj = max | — [a n U™+i + (1 - a n )U^] , 

L Pn J 


n = N — 1, ■ ■ ■ , 


1 , 0 . 


(27) 


Now we consider the relation of BTM and explicit difference scheme for American option. 
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Lemma 3 If lnu = Ax = o n \JAt n , we have 

Here 9 n are coefficients of BTM defined by ([ 6 ]). 
The proof is easy. 


Contrasting (101 and (27), BTM is equivalent to a special explicit difference scheme (27) 
in the sense of neglecting 0(Ax 3 ). 


Now we show the conditions for American (put) option price be monotonic. 


Theorem 4 Assume that 0 < a < 1 and 



(i) If tpj = {S^ Ax — E) + (call), then Uf < U(f +1 

(ii) If Lpj = (E — Soe^ Ax ) + (put), then Uf > f7 " +1 


n Qn 



< 1 . 


and 0 < Uf 1 < e jAx+c . 
and 0 < Uf < E. 


Proof From the assumption we have 0 < a n < 1. 

(i) Uj 1 = ipj = (Soei Ax — E) + < (S'oe^ +1 ^ Ax — E) + = tfij+i = Uf +1 . Now assume that 

Uk+i < V k+i_ Then we have 

Uf = max |i- {(1 - a)Uj +1 + a [a k U*+) + (1 - a k )U^] } , j 

< max {1 {(i - a)U^l + a [a k U$ + (1 - a k )U* +1 ] ) , ^ +1 J = U? +1 . 

(ii) is proved in the same way as (i). (QED) 


Aa( r _ < 

a 2 n y™ 2 


^ < 1 , r(t)/a 2 (t) is increasing 


Theorem 5 Assume that 0 < a < 1, 

and q(t)/a 2 (t) decreasing on t. Then prices Uf of American put option given by (25) and 
(26) with ipj = (E — So^ Ax ) + are decreasing on t, that is, 

U™>U™+\ jeZ, n = N, 1 , 0 . 


Proof When n = N — 1 , from (261 we have Uj > ipj = U A , j € Z . Assume that 
Uj +1 > Uj +2 , j £ Z . From the assumption and lemma [I] (i) we have p k < p k +i and thus 

Uf = max | {(1 - a)U* +1 + a [a k U*+) + (1 - a k )U^+l] } , | 

l ~ {(1 - «)U* +2 + a [a k U*+? + (1 - a k )U^ 2 ] } , ^ j . 


> max 


From a k = \ + and the assumption, we have a k < Ofc+i • By theorem^ 

(ii), we have Uj+ 2 > U^ 2 and thus lemma [ 2 ] with a k < a k + 1 gives us 


a k+ iU^+ 2 + (1 - a k+1 )U£? < a k U*+( + (1 - a k )Upf. 


rk-\-2 


k +2 


rfc+2 
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Therefore we have 

Uf > max (— {(1 - a)Uj +2 + a [a k+1 U*+? + (1 - a k+1 )U^+ 2 ] } , ^} = U* +1 . 

I Pk+l J 

(QED) 

Remark 6. Theorem [5] strongly represents the effect of time dependent coefficients. 
Here the main tools are lemma[l](i) and lemma[2j The conditions of theorem[5]are essential. 
If r(t)/a 2 (t) is not increasing, then the price of American put option by explicit difference 
scheme might not be decreasing on t as in remark 4. 

Remark 7. Only using the analogs of lemma |T] and lemma [5] it seems difficult to prove 
that American call option’s price is decreasing on t. 

Now we show the existence of approximated optimal exercise boundary. 

Theorem 6 Under the assumptions of theorem [dj for any 0 < n < N ~ l , there exists 
j n G Z such that 

j < jn =>• U™ = j = j n + 1 => Uf > <Pj\ j >jn + 2=> U.J > tfij. (28) 

jo < 31 < ••• < jN-1. (29) 


Proof Note that t/A = (E — 5'oe ; ' A:I ’) + = is decreasing on j G Z. Let 

fci = max{j G Z;E — S' 0 e' jAx > 0}. 

Then if j < k\ — 1 then j — 1, j, j + 1 < k\ and 

(Pi -1 =E- S 0 e jAx ~ Ax , (fij = E — S 0 e jAx , ip j+1 = E — S 0 e jAx+Ax > 0. 
Let u = e Ax , d = e~ Ax and ifj = - S 0 vP ■ then 

’ ' J PN—1 U fiKT— 1 


pN-1 


U. 


N-l 


= max 


1 —— [(1 ^ a)(pj + a(a N -itpj + i + (1 - ajv-i)<Pj-i)] > <Pj \ 
l Pn-i J 


= ma ifij}. 


Then we have 


E 


lim if)j = - < E = lim ipj, (j < k\ — 1). 

j^y-oo Pn-i j-y-oo 


(30) 


First, we consider the case that ipj < ipj (Vj < k\ — 1). If j < k-[ — 1 then U- = ipj, 
and if j = k± + 1 then ipj = ifij+i = 0, <Pj—i > 0 and thus we have 

a(l — a^-i) 


Uf- 1 = max 


Pn-i 


If j > k-[ + 2 . then from (261 we have U A 1 > <pj . So if t/A 1 > tpj k , then we define 


<Pj- 1, 0 > > 0 = ipj. 


N -1 


jjv_i = ki — 1 ; and if 1 = ipj h , then we define jjv-i = k\ . 
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Next, we consider the case that 3 j (j < k\ — 1) : ipj > Pj ■ We define 
j n —i = niax-jj < fci - 1 : ipj < Pj}. 


Then we have 


j — j N— l ^ Pj — [ Pj ^ Uj — tpj , 

j = j N — 1 + 1 =>“ 'Pj > '-Pj ^ Uj = pj > tfj , 

3 > Tv-i + 2 ^ ipj. 

Thus we proved the existence of jjv-i < &i. 

Now we assume that when n = k + 1 there exists jk+i such that 


jk+l < jfc+2 < 
j S ife+i = 
j = jfc+i + 1 = 
j > jfc+i + 2 = 


< JN-1, 

Pj, 


U j +1 = 


U * +1 > Pj, 


Uj +1 


If j < jk+i - 1 , 

the above, let ip 3 

U t 


then j 

E 


+ 1, j; 


> (flj. 

k +1 


= it - 


j - 1 < jk+l and thus t/f +i = < p t (i = j - 1. 

l- a+a [a k u+(l-a k )d] _ Then by Q we have 


(31) 

3, 3 + !)■ As 


= max« 


; { T- [(! — Oi)pj + a(a k p j+ 1 + (1 - a k )<pj -\)\, 

= max{^j, <pj}. 

Note that ipj < <pj for sufficiently large j £ Z. In the case that ipj < Pj (Vj < j k +i — 1), we 
have Uj = <pj for all j < j k +i — 1 • From theorem |s] and the inductive assumption (31) we 


u*>u 0 


k +1 


> Pj- 


have the fact that j = j k + 1 + 1 => Uj > Uj > <pj; j > j k + i + 2 

Therefore if E /£ +1 > p jk+1 , then let j k = j k +1 - 1 . If U!f k+1 = p jk+1 , then let j k = j k +1 . 
In the case that 3 j (j < j k +\ — 1) : ipj > Pj, we define j k = max{j < j k +i — 1 : ipj < tp^-} • 
Then 

i ^ ifc br = Pj, 


3= jk + l - l j 
j + jk + 2 =+ L j 


Uj = ipj > Pj > 
Uj > Pj- 


Thus we proved the existence of j k < jk+i- (QED) 

Remark 8. If Ax is enough small, then j k £ [jk+i — 1, jk+i] ■ 

Now we estimate the optimal exercise boundary near the maturity. 

In the first part of the proof of theorem [ 6 j we proved the existence of jfv_i 
approximated optimal exercise boundary near the maturity. If k\ is the one defined in (30) 
and Sq = e c , then k\ = max-jj £ Z;E — e^ Ax+c > 0} and for j < k\ — 1 , we have 
Pj = E — e i Ax+c and let 


the 


iPj = 


1 


Pat-i 


[(1 - a)pj + a(a N _iPj + i + (1 - a N ^)pj_ i)]. 
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In the case that ipj < tpj (Vj < k\ — 1) we know jlv-i = fci — 1 or k\ . Then we have 

E — gJw-iAx+c > o, — e (Jjv-i+2)Ax+c < q anc j ^us we j mve 

In if — 2Ax < jJv-iAx + c < In If. (32) 

In the case that 3 j ( j < fci — 1) : ipj > <Pj, by theorem [6j we have 
j n —i = max{j < ki - 1 : ipj - ipj < 0}. 

By using the definition of a n and Taylor expansion, we have 


a n e Ax + (1 - a n )e Ax = 1 + gn Ax 2 + 0(Ax 4 ). 

< 


Then for j < k\ — 1 we have 


Vh - = 

(1 - a)(E - e jAx+1 ) + a [a N -i(E - e 0' +1 ) Ax+1 ) + (1 - ojv-i)(-B - e (j- 1 ) A ^+ 1 )] 


-(£ - e JAx+1 ) = 


1 (jj^_^Ati^—i 


Pn-i 


Ax 2 


PN -1 

At 2 

fev_ ie ^ +c - + 0(Ax 4 ) 

CT Af-l 


Note that E > e^ Ax+c for j < k\ — 1. If q ^~i < rjv -i , then r^-iE > q^~ 1 e : > Ax+c and 
therefore if Ax is enough small, then we have ipj < ipj (Vj < k\ — l) . Thus in our case, since 
3j ( j < fci — 1) : ipj > ipj, we must have qN-\ > rjv-i • Then for sufficiently small Ax , we 
have jjv-i = max{j < k\ — 1 : qN-i^ Ax+c < tat-i-B} and therefore jv-i Ax+c < In ^A=^E. 
= j n —i + 1 , then qN~iej Ax+c > tn-iE and thus {jn-x + l)Ax + c > In E . So 


If 3 = 3N 

we have 


In 


nv-i t-, a . A , ,, rv_i 

-If — Ax < jat-iAx + c < In-If. 


qN-i 


qN-i 


Thus combining this inequality with (321, we have the following theorem which provides 


an estimate of the approximated optimal exercise boundary near the maturity. 


Theorem 7 In min 


in ( E > ~ 2Ax - 3 n-\Ax + c < In min (e, • 


For fixed Ax , the approximated optimal exercise boundary x = pAx{t) is defined as 
follows: 


PAx(t) = 7 ^—^A_(j n+1 Ax + c) + — j-(j n Ax + c), t G = 0, ■ • ■ ,N - 2. 


^n +1 


^n +1 


Corollary (i) PAx{tN-i) e [inmin (if, ^-^If) - 2Aar, lnmin (if, j-.E 1 ) 
(ii) PAx{t) is increasing ont. 
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5 Convergence of the Explicit Difference Scheme and BTM for 
American Put Option. 


In this section we will prove that the explicit difference scheme (25) and (26 1 for American 


put option converges to the viscosity solution to the variational inequality (181 and using it 


prove the monotonic properties of the price of American put option and its optimal exercise 
boundary. 

We denote by l°°(Z) the Banach space of all bounded two sided sequences with sup 
norm. In l°°(Z), we define (■ Uj ) < (Vj) <t=> Uj < Vj,\/j £ Z. 

For fixed every n, the two sided sequence U n = (••• , Uj, ■ ■ ■ )'jL_ 00 of American put 
option’s prices Uj , Vj £ Z given b y ([25] ) and (261 is bounded from theorem [4j 

If we denote the right side of (26) by (F„U n+1 ).,- , then F„ defines an operator 


U" := F n U " +1 = {(F„U" + 1 )j}^ : _ 0O 


(33) 


sending the sequence U " +1 of t n+ \ -time prices to the sequence U" of t n -time prices. The 
operator F n depends not only on n and Ax but also on t n and A t n . 


Lemma 4 If 0 < a < 1, 


Ax 


[r n -Qn- 


< 1, then F n is increasing, that is, 


U < V, U, V G r(Z) => F„U < F„V. 

(Proof) Noting that from the assumption we have 1 — a > 0 and 0 < a n < 1, the 
required result easily comes from ([26]) . (QED) 


Lemma 5 If U £ l°°(Z ) and K = (• • • , K, K, K ■ ■ ■) with I\ > 0, then 

F n (U + K) < F„U + K. 

(Proof) Since p n > 1 , we have 
F„(U + K) = 

= (max | — [(1 - a)(Uj + K ) + a (a„(U j+1 + K ) + (1 - a n )(Uj-i + K))\, ipj 

\ l Pn 

< (max | — [(1 — a)Uj + a (■ a n Uj + \ + (1 — a n )Uj-i)] , ipj \ + max < —, 0 
V l Pn J L Pn 

<F„U + K. 


j=-oo 


j = -oo 


(QED) 


Define the extension function u^ x {x, t ) as follows: 

x £ [(j - l/2)Ax + c, (j + 1 /2) Ax + c),t £ [t n ,t n+1 ) => uax(x, t) := Uj, (34) 
x £ [{j - 1/2)Ax + c, {j+ 1/2)Ax + c),t £ [t N ,T) =>• u Ax (x, t) := Uf. (35) 

Remark 9. Here max is piecewise continuous function. The following discussion is also 
true if we define uai as a continuous function interpolating the data set {jAx + c, [/") . 
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From Theorem [4j we have 

0 < ua x {x, t) < E. (36) 

Therefore, for every fixed t, if we denote t) := {u/\x(x, t) : x € R}, then we have 

UAx(*, t) G l°°(Z). 

When t G [t n , t n+ 1 ), n = 0, • • • , N — 1, if we define 


At = At(t, Ax) := ^—%-Ai n+ i + * n+1 * At n , 


t"n +1 In 

then t + At £ \t n +\, L 1 + 2 ) and thus we have 

UA*(», t) = F n UAx{*, t + At). 


tn +1 t n 


(37) 


(38) 


A t„ 


. . ^ + . ** +1 t • 1 is a convex linear combination of -°^r~ and 1. Thus if 

tn+l—tn CT“ tn + l—t n <T“ 


<j(t) is continuous, then this ratio converges to 1 when Ax —> 0 . Like this, we proved the 
following lemma: 

Lemma 6 Assume that cr(t) is continuous. Then limAz-^o At ^ Ax 'i = 1 . 

In order to prove the convergence, we recall the notion of viscosity solutions. 


Let U S'C([0, T\ x R) (LS , C([0, T] x R)) be the space of all upper (lower) semi-continuous 
functions defined on [0,T] x R. If u G USC([0,T] x R) (LSC([0,T\ x R)) satisfies the 
following two conditions, then u is called the viscosity subsolution (supersolution) of the 
variational inequality ( |18| ): 

(i) u(x, T) < (>Mx), 

(ii) If $ G C 1 , 2 ([0,T] x R) and u — $ attains its local maximum (minimum) at (x,t) G 
[0, T] x R , we have 


mm < —-- 


5$ a(t)' 2 


dt 


2 dx 2 


. a(t) 2 \ 9$ 

- ( r(t) - q(t) - — j —+ r(t)u, u - 


<(>) 0 


(x,t) 


u G (7([0,T] x R) is called the viscosity solution of the variational inequality (18) if it is 
both viscosity subsolution and viscosity supersolution of (18). 


Lemma 7 (Comparison lemma) [7 If u and v are viscosity subsolution and supersolution 
of (18) and \u(x,t)\,\v(x,t)\ < E, thenu<v. 


Theorem 8 |7j/ Ifr(t), q(t) are continuous on [0, T], then the problem fl8| has a unique 
viscosity solution. Furtheremore, there exists an optimal exercise boundary pit) : [0, T] —> R 
(continuous function) such that if x < p(t), then u(x,t) = <p(x) ; if x > p(t), then u(x,t) > 
p{x) and in this region u{x,t) is a classical solution to the equation 


du _ (j{t) 2 d 2 u 
dt 2 dx 2 


(t) ~ v(t) 



du 

dx 


+ r(t)u = 0 . 
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Remark 10. It is easy to show p(T) = ln£Jmin[l,r(T)/g(T)], using the way of 10 . 


Theorem 9 Suppose that u{x,t) is the viscosity solution of (18 1 . Assume that 0 < a < 1, 
r(t)/a 2 (t) is increasing and q{t)/cr 2 (t) decreasing on t. Then we have 

(i) UAx{x,t) converges to u(x,t) as Ax —» 0. 

(ii) pAx{t) converges to p{t) as Ax —> 0. 


Proof Suppose that it(x, t) is the viscosity solution of (18 1 and denote 

u*(x,t) = lim sup UAx(y, s), 

Ate—>-0, (y,s)—^(x,t) 

u*(x,t) = lim miuAx(y,s) 

Ate—>-0, (y,s)—>(x,t) 


(39) 


From (36), u* and it, are well defined and we have 0 < u*(x,t) < u*(x,t) < E . Obviously 
u* £ U 5C([0, T\ x R) and it* £ LSC{[ 0, T\ x R). If we prove that u* and it* are subsolution 
and supersolution of (18), respectively, then from lemma [7J we have u* < u * and thus 
u* = it, = it(x, t) becomes a viscosity solution of (jTsj), and therefore we have the convergence 
of the approximated solution UAx{x,t) . 


We will prove that u* is a subsolution of (18). (The fact that it* is a supersolution is 
similarly proved.) From the definition (35), we can easily know that 


u*(x,T) = tp{ x) = (£-e x )+. 

Suppose that £ C 1,2 ([0, T\ x R) and u*—(j) attains a local maximum at (xo, t 0 ) £ [0, T) x R. 
We might as well assume that (it* — <t>)(xo,to) = 0 and (xo,fo) is a strict local maximum 
on B r = {(rr, t) : t 0 < t < t 0 + r, |x — Xo| < r}, r > 0 . Let $ = (j) — e,e > 0, then u* — $ 
attains a strict local maximum at (xo,to) and 

(u*-$)(x 0) to)>0. (40) 

From the definition of it*, there exists a sequence UAx k {Vk, s k) such that Ax* —> 0,y k —t 
x 0 , s k -T t 0 and 


lim UAx k (yk,Sk) = u*(x 0 ,t 0 ). (41) 

k—foo 

If we denote the global maximum point of UAx k — $ on B r by (y*,, Sk), then there exists a 
subsequence UAx k . {ilk, Sk) such that 

Ax ki ->■ 0, ( : yki,S ki ) -> (x 0 ,t 0 ), («A It . - ^){Vki, S ki ) -t ( u * - $)(x 0 ,t 0 ) as h ->• oo. (42) 
Indeed, suppose (y*^, S ki ) —» (y,s), then from ( |4l| ) we have 

(it* - <F)(x 0 ,fo) = lim (uAx k . -®)(yki,Ski) < lim {uAx k . - §){y ki , S ki ) < (u* - <F)(y, s). 

ki~t oo 1 ki~> oo 1 


Therefore we have (y, s) = (x’o, t 0 ), since (xo, to) is a strict local maximum of (it* — $). Thus 
for sufficiently large k *, At = At(Ax ki ) defined by (37) is small enough and if {x,S ki + 
Af(AxfcJ) £ B r , then we have 


{uAx k . ~ 4>)(x, + At(Axki)) < (UAx ki ~ ®)(yki, S^), 
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that is 


UAx k . (x,S ki + At(Ax ki )) < $(x,S ki + At(Ax ki )) + (u A x k . ~ $)(y ki ,S k .). (43) 


From (40) and ([42]) , we have 

{uAx k . - ®){yki,S ki ) > 0 (for sufficiently large hi). 


(44) 


aAxi 


For every ki, define t n and A t n = a2 ^ t ) as in (19 1 with Ax ki . Select t ni and j ki = j such 

that y k . G [(j - l/2)Ax ki + c, (j + l/2)Ax fci + c) , S ki G [t ni ,t ni + i), and simply denote 
tm = t n and j = j k . ■ Then from ( |38[ > , lemma [4] with ( 431 and lemma [5] we have 

uax k {y ki ,S k J = U J* = (F n U n+1 )j = [F n u Atkt (., S ki + At(Ax k MykJ 

< {F n [&{»,S ki + At) + (uAx k . - ®)(yki,S ki )]}{y ki ) 

< {F n [&(»,S ki + At)]}{y ki ) + (uAx ki ~ ®)(y ki ,S ki ). 

Thus we have 

$(y ki ,S ki ) - {F n [^(*,S ki + At)]}(y ki ) < 0. 


Therefore using (26) we have 


®(y ki ,S k J - max 


1 


1 -{- v n At n 


1 ~ + At ( Ax ' fc .))+ 


+ °A^ n [ a n$(y ki + A x kil S ki + At) + (1 - a n )$(y ki - A x ki ,S ki + At)] 
This inequality is equivalent to the following. 


> m < o. 


Ai„ 


1 -(- v n A t n 


^iVkuSki) - $(y ki ,S ki + At(Ax ki )) 
A t n 


°n ^{Vkj + , 5 fci + At) - 2<F(y fc ., Ski + At) + ®(ykj - Ax ki , 5 fci + At) 

2 2Ax k , 


r n - 


4>(g fci + A x ki ,S ki + At) - $(y ki - Ax ki ,S ki + At) 
2Ax ki 


+r n <b{y ki , S ki ) , $(y ki ,S ki )~ tfj j < 0. 


Noting that 

I 


l+r^At,, > °> We haVe 


. n / ®(ykj,S ki ) - $(y ki ,S ki + At(Ax ki )) 

A t n 

T n ®(yki + Ax ki , S ki + At) - 2 ${y ki , S ki + At) + 4>(y fci - Ax ki , 5 fci + At) 
2 2Ax ki 

&n\ $ ( y ki + A x ki ,S ki + At) - <F(y fci - A x ki ,S ki + At) 

2Ax k 


- 1 r n -q n - — 


+r n ®{y ki , S ki ),®{y ki , Sk.) - v’jj < 0. 
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This inequality is equivalent to the following. 

$(y ki ,S ki ) - $(y ki ,S ki + At(Ax ki )) 
At.(Ax ki ) 


At(Ax ki 
A t n 


®(ykj + Ax ki , S ki + At) - 2$(y ki , S ki + At) + $(y ki - Ax ki , S ki + At) 
2 2Ax kt 

cr 2 \ $(y ki + Ax ki ,S ki + At) - <&{y ki - Ax ki ,S ki + At) 


~ f n ~ Qn - 


2Ax k 


+r n $(y ki ,S ki ),$(y ki ,S ki ) - <Pj} < 0. 
oo, then Ax ki —> 0. From lemma j^J we have At ( Axk ^ 


Let ki 


A 


1. Thus we have 


mm < —-- 


9$ 

dt 


2 (t) 9 2 $ 
2 dx 2 


cr 2 (t)\ 94> 


- ( r(t) - q(t) -— ) — + r(t)$, $ - <p 


< 0. 


(Here we considered ( y ki , S ki ) —► (xq, to) and Pj k 


(®o> *o) 

<p(xo)- Let e —► 0 , then we have 


d<t> 

mm <, — — - 


er 2 (t) d 2 4> 
2 9a; 2 


- »*(t) - q(t) - 


<t 2 (t)\ dc/) 
dx 


+ r{t)<t>, (/)-ip 


< 0 . 


(*o, to) 


Since u*(x o, to) = </>(a:o, to), u* is a subsolution of (18 1 . Thus we proved (i). 

Now we will prove pAx{t) converges to p(t) as Ax —> 0. The main idea is from 12 
First, from Corollary of theorem [7] and Remark 9, we have 


lim pAx{tN-i) = In min 
Ax —>-0 


. ’q(T) J 


= P(T). 


Now fix to £ [0,T) and suppose x < limAx->osuppAx(to)- Then there exists a sequence 
Ax k such that Ax k — > 0 and limfc->oo PAx k (to) > x. Denote by {t^" 1 } the time partition 
corresponding to Ax k and let to £ [t^_i, Then tn k J — > to as k —> oo. Since pax is 

increasing, we have 

PAx k {tn'J ) > PAx k (to) ■ 

Select x k such that lim Aa;*, = x. For sufficiently large k, x k < pAx k {trio) = jn 0 Ax k + c 
and thus we have UAx k ( x k itrio) = <f{x k ) (since x k is in the exercise region). Thus we have 
Ax k —»• 0 => u{x, to) = ip(x) => x < p(to), so we have 


lim sup pax (to) < p(to)- 


Now we will prove liminf pAxiko) > p(to). Assume that there exists e > 0 such that 
lim inf pAx(to) < p(t o) — 2e . From the fact that p(t ) is continuous, there exists S > 0 such 
that liminf pax (to) < p(t) — 2e,t £ [to — 6 , to + 5] ■ Therefore there exists a sequence Aa;*, —> 0 
such that pAx k (to) < p(t) — 2e, t £ [to — <5, to + £]• Now let p := min{p(t) : t £ [t 0 — 5, to + <5]} 
and Q = {{x,t) : t £ (to — S,to],x £ (p — 2e, p — e)}. Then since pAx k {t) is increasing, we 
have pAx k (t) < PAx k (to) <p — 2e < x < p ~ e for (a:, t) £ Q. Therefore we have 

PAx k {t ) < x < p(t) - e, (a;,t) £ Q. (45) 
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From pAx k (t) < x, let UAx k {x,t) = £/", then we have 

U ? +1 ~ U™ a 2 (t n ) - 2U ? +1 + U ?-i 


At 


r(t n ) - q(t n ) - 


cr 2 (t„) 


Ax 2 

Tjn+l _ rrn+1 

U j + l U j -1 

2Ax 


+ r{t n )W = 0. 


Letting Air*, —t 0 , then UAx k {x,t) (= [/") converges to the viscosity solution u(x,t) and 
from the regularity of viscosity solution (Theorem [ 8 ]), we have 


du a(t) 2 d 2 u 


m + 


2 dx 2 


+ r(t) - q(t ) - 


er(t) 2 \ <9u 


— - r (t)u = 0, (x,t) G Q. 


(46) 


On the other hand, from (451, we have x < p(t ) — e , thus x is in the exercise region of (18) 
and u(x, t) = E — e x . So we have 


du 

m 


r(t) 2 8 2 U 

2 dx 2 


+ r{t) - q(t) - 


<t( t) 2 \ du 
dx 


- r{t)u = q{t)e x — r{t)E , (x, t ) G Q. 


This contradicts to (46). Thus we proved 

liminf/9Ax(to) > p(t 0 ) > limsup pAx(t 0 ). 


So we have limpAx(io) = p(^o)-(QED) 


Corollary 1 (Monotonicity of American put option’s price and optimal exercise boundary) 
Suppose that u(x,t) is the viscosity solution of (18). Assume that r(t)/a 2 (t) is increasing 
and q(t)/a 2 (t) is decreasing on t. Then we have 

(i) u(x,t) is decreasing on x and t. 

(ii) p(t) is increasing on t. 


(Proof) (i) comes from Theorem [4] (ii), Theorem [5] and Theorem [9] (i). 

(ii) comes from (ii) of Corollary of Theorem [T] and Theorem IU(ii)- (QED) 


Lemma 8 1A] Let fl C R m and f n (x i, • • • ,x m ) be pointwise convergent to a continuous 

function f{x i,--- , x m ). Assume that f n and f are monotone on fb Then f n uniformly 
converges to f on any compact subset of Q. 

Theorem 10 When Ax —> 0, then UAx{x,t) uniformly converges to u(x,t) on any bounded 
subdomain of[0,T] x R and pAx(t) uniformly converges to p(t). 

(Proof) From the result of Theorem [To} u(x, t) and pit) are both monotone. Thus from 
lemma |8j we have the desired results. (QED) 

Remark 11. Consider the convergence of the binomial tree methods. As shown in 
section 3, in the lattice Q c = {( Xj , t n ) : Xj = j Ax + c, 0 < n < N — 1, j G Z}, the 
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explicit difference scheme (27) with a = <r 2 At n /Ax 2 = 1 coincides with BTM in the sense 
of neglecting 0(Ax 3 ). 

Let V " be the prices by the explicit difference scheme (27) and VJ 1 the BTM prices. 
Then we have 

U J* = max ( — [a n U^ + (1 - a n )U , (E - S 0 e jAa: ) + } , 

L Pn J 

Vj 1 = max { — [O n V£tf + (1 - 0 n )V. "ft 1 ] , (E - S 0 e jA *) + l . 

I Pn J 

Let u = e Ax . Then we have 

\vp -u?\<~ [I e n v?£ a n w $| +1(1 - 0 n )vp? - (1 - a„)c/;_+ 1 |] 


Pn 


1 


< - [o n \vjjj 1 - up+j\ + K - a n \up+j + (i - o n )\vp? - u;+j\ + \e n - a n )\up+j] 

Pn 

< — [||V” +1 - U" +1 ||^ (z) + 20(Aa; 3 )||U n+1 11(00(2)] . 

Pn 

Here we considered lemma |3j 9 n = a n + 0(Ax 3 ). Note that ||U n+1 ||(oo( 2 ) < E, then we 
have 

\v?-u?\ < — [||V" +1 — U Tl+1 ||;oo ( 2 ) + 2E • 0(Ax 3 )] . 

Pn 

Since r^/cr 2 is increasing, we have 

exp(-r„Af n ) = exp(—r„ A.T 2 /cr 2 ) < exp(-r 0 A:c 2 /cTo) =: A, 
and using p ” 1 = (1 + r n At n ) = exp(— r n At n ) + 0(At 2 ) , then we have 
|V7 - t/J 1 I < ri||V” +1 - U” +1 11(00(2) + 2 AE ■ 0( Ax 3 ). 


Therefore we have 


|| V” — U” |(o.( 2 ) < 

H||V ” +1 - 

< 

A 2 v n+2 - 

< 

A iV- n || V iV 

Here noting that || — 

^11(00(2) = 


D (z) 


Tn+2| 


( 00 ( 2 ) + (A + A 2 )2E ■ 0(Aa; 3 ) 


V N - U" 11 ( 00 ( 2 ) + (A + A 2 + --- + A N ~ n )2E ■ 0(Ax 3 ). 




A 


IIV" - Ul(oo ( 2 ) < Izrj 2E ■ 0(Ax 3 ) = 0(Ax). 

Thus from theorem [9j when Aa; —> 0, the BTM prices Vj 1 converges to the viscosity solution 
to the variational inequality (18). 
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